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Abstract: In this letter, we present the NNLL-NNLO transverse momentum Higgs dis¬ 
tribution arising from gluon fusion. In the regime p± <C mn we include the resummation 
of the large logs at next to next-to leading order and then match on to the a% fixed order 
result near p± ~ my- By utilizing the rapidity renormalization group (RRG) we are able to 
smoothly match between the resummed, small p± regime and the fixed order regime. We 
give a detailed discussion of the scale dependence of the result including an analysis of the 
rapidity scale dependence. Our central value differs from previous results, in the transition 
region as well as the tail, by an amount which is outside the error band. This difference is 
due to the fact that the RRG profile allows us to smoothly turn off the resummation. 
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1 Introduction 

We are transitioning into an era of precisions Higgs physics. Presently there is sufficient 
data to study various decay modes of the Higgs [1], and soon there will be sufficient data to 
study the the Higgs differnential cross section, which can be used to better understand the 
underlying production mechanism and to search for new physics [2], While the new physics 
would be most prominent at large values of p±, the predominance of the events will be in 
the lower range. Furthermore, to isolate the Higgs’ decays from backgrounds, events are 
binned according to their highest-transverse momentum jet. The 0-jet bin, corresponding to 
no central jets above a certain transverse-momentum threshold, plays a critical role in the 
current Higgs analysis 1 . Thus there is increased motivation for making precise predictions for 
the distribution when p± <C mh since back to back central jet production is power suppressed 

[15]- 

1 Resummed predictions for the 0-jet bin can be found in [4-8]. 
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The small p± region of parameter space is polluted by large logarithms which must be 
resummed in order to retain systematic control of the theoretical errors. Resummations have 
been previously been discussed within an SCET [3] framework [[17],[20],[21],[18],[22],[23], [25], 
[27], [26]] as well as in the CSS resummation formalism [[28],[29],[30],[31],[33],[32]]. While 
formally most of the these results agree at a given order in the resummation procedure, the 
central values as well as the predicted errors between different resummation techniques will 
differ - ’. Another crucial difference between various theoretical predictions is how the transition 
between the fixed order result at large p± and the resummed result is handled. In this paper 
we present a result for the resummed cross section at NNLL + NNLO within the confines 
of SCET and the rapidity renormalization group (RRG) [14, 15]. Our motivation for this 
analysis is two-fold. By working within the RRG formalism we are able to consistently turn 
off the resununation in the tail region. It has been shown in the context of B —> X s + 7 
[48], thrust [49] as well as the Higgs jet veto calculation [44], that if one does not turn off the 
resummation then one can over estimate the cross section, in the region where fixed order 
perturbation theory should suffice, by an amount which goes beyond the canonical error band 
in the fixed order result. This overshoot happens despite the fact that the resummed terms 
are formally sub-leading in the expansion. The reason for this overshoot has been shown 
[44, 48, 49] to be due to the fact that there are cancellations between the singular and non¬ 
singular terms in the tail region and that this cancellation will occur only if the proper scale is 
chosen in the logarithms. This will be born out in our analysis as well. By using the RRG we 
have a natural way to interpolate between the resummed and fixed order regions. This allows 
us to present the full spectrum such that in the large p± region our results match onto the 
NNLO result, as previously discussed within the context of the jet veto cross section [7, 8 ]. 
In this sense our paper fills a gap in the literature. As we will discuss below, the RRG plays 
a important role in the theory error determination. We also discuss in detail how our error 
estimates and full spectrum differ from those of previous work. 

2 Systematics 

In this section, for the sake of completeness, we review the well known systematics of the 
calculation. We work within the confines of the large top mass approximation. This approx¬ 
imation is known to work extremely well - much better then one would naively expected - 
especially away from the tail of the spectrum. We will be not be considering p± large enough 
for these corrections to be relevant to the error budget. A full discussion of such errors as 
well as others will be discussed at the end of the paper. 

At the scale rrq, full QCD is matched onto the large top mass effective theory at order a 
after which we match onto SCET at the scale rrih- In SCET we will work to leading order in 
a systematic expansion in p±/rrih, and the errors due to non-perturbative corrections, which 
are suppressed by A/pj_, are included in the final error analysis. 

2 For a detailed comparison of resummation methods in the context of e + e~ event shapes see [57], and for 
threshold resummations see [34, 35]. 
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As is well known, logs in the perturbative expansion - in impact parameter space - 
exponentiate and allow us to organize the series as follows 


(7 ~ 


+ a s H(a s In 


+ .... 


( 2 . 1 ) 


The standard terminology is such that keeping only F corresponds to “LL” leading (double) 
log, while if G is retained we would have NLL etc. As p± gets larger the EFT breaks down 
and the fixed order calculation becomes the relevant quantity. We will utilize the order a 2 
(NNLO) result [9], [10] 3 in this large p± region. We do so despite the fact in the small p± 
region we will be only keeping terms of order a. This is not inconsistent as the error in 
the two disparate regions are distinct and uncorrelated. In the resummed region we will be 
working at NNLL in the perturbative expansion and at leading order in the power expansion 
V^l m h- By utilizing the RRG scale we will smoothly turn off the resummation and match 
onto the full NNLO fixed order calculation. 


3 Factorization in SCET and anomalous dimension 


A factorization theorem for Higgs production at small p± within SCET and the RRG for¬ 
malism was developed in [15]. The starting point in the large top mass effective theory is the 
gluon fusion operator 

U(x) = C t h ^-Tr[G^(x)G fiU (x)\. (3.1) 

The matching coefficient for this operator ( Ct ) is known to two loops [12]. 

At the scale 4 we match onto SCETII. This theory is the version of SCET in which 
the collinear modes and the soft modes have the same invariant mass, which distinguishes 
it from SCETI where there exists a hierarchy in masses between these two modes. Due to 
this equality in virtualities a new set of divergences - rapidity divergences - arise, which are 
not regulated by dimensional regularization. The rapidity regulator introduces a new scale 
which acts as a boundary between the collinear and soft modes as discussed in [14, 15]. The 
reader may consult ([15]) for the details of the formalism. At leading order in A = p±/Mh the 
differential cross section for higgs production at low transverse momentum may be written as 


da 

dp]_dy 


C? 


8 v 2 S(N? - 1) 
4(2vr) r 


d xe 


^^(2ir)5 + (pl-ml)5^y - &(p± ~ \Ph ±\ 2 ) 

lx ' Ph H(m h )f^ g/p { 0, x + , x±)f ±g/Pf , l/ (x~, 0, x_l)cS(0, 0, x±) 


(3.2) 


3 As is common with fixed order calculations, the authors of [10] use the term NLO, since the leading order 
result is a delta function in p±. 

4 We ignore the running between the top mass scale and mu as these logs will be sub-leading in our power 
counting. 
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S is the soft function defined as 


5(0, 0, x±) 


1 

(27r)2(JV c 2-l) 


{0|Sr(x)S? J (x)S;“(0)Sr(0)|0), 


(3.3) 


defined in terms of the light-like Wilson lines 


S n (x) = Pexp ( ig n ■ A s (nX + x)d\ 
\ J — OO 


(3.4) 


f a/S is the transverse momentum distribution function (TMPDF) which is matched onto the 
PDFs at the scale p± 2> A. 


/f o/p (0,x+,x ± ) = oToTvi (Pn\[B^(x + ,x±)B^{0)]\p n ), 


’“'^“2(27r) 

f±g/p( X ~,°,X±) = 9(27r)3 (P ft l ( X ~» f (0)] |Pfi) 


(3.5) 


defined in terms of 


= -Tr \T a \Wkx)iD» ± W n (x) 


9 


W n is a collinear Wilson line in the fundamental representation defined in x-space by 


W n (x) = Pexp ^ ig J n ■ A n (n\)d\ S J . 


(3.6) 


(3.7) 


n/n = (1,0,0,±1) are the null vectors which correspond to the directions of the incoming 
protons. In impact parameter space, we can write the functions S and f a P in terms of their 
inverse Fourier transform. 

s(b) = J WF e<lF;±s{p:±) 

/r(M = I ^ e&Pl/r(p± ’ z) - 

The cross section is given by 

ijfi V = ^ 2l ' f 2^mi0T) / < 3 - 10 ) 

H(Mh ) is the hard coefficient which depends on the scale p = Mh and and we define H(Mh) = 
8M?|C' S | 2 , where the coefficient C s is known to two loops [13]. The soft function as well as the 
TMPDFs depend upon both p, the usual scale introduced within the context of dimensional 
regularization which factorizes the hard modes from the soft and collinear modes, as well as 
v, the rapidity factorization scale which distinguishes between the soft and collinear modes. 


(3.8) 

(3.9) 
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4 Renormalization and Resummation 


Each of the pieces of the factorization theorem H, S and /, have natural scales with which 
they are associated. H is independent of v and its natural p scale is nih- The soft and 
collinear functions have natural scales (p = p±,v = p ±) and (p = p±,v = rrih) respectively. 
When these objects are evaluated at these scales, they will be devoid of large logarithms. 
However, given the natural distribution of scales we can see that it is not possible to choose 
a p and v such that all of the individual pieces sit at their respective natural scales. Thus 
it is expedient to choose a (p, v) value such that the maximum number of pieces sit at their 
natural scale, and then use the RG and RRG to sum logs for the remaining pieces. 

The individual functions satisfy the following RGE’s. 

= tfH(p) 

Ki ~dL' F ( b, p ,u ) = 7 £ F ( 6 ’t 4 - 1 ) 

(4.2) 

where Ki = (p, u) and F E {5, The anomalous dimensions satisfy the relations 

7 ? + 7 ? + 2 7 1 = 0 

7 v + H = 0. (4.3) 

All the renormalization group equations are diagonal in impact parameter space and hence 
are straightforward to solve. The scales for RG and RRG operations commute, and as a 
consequence of this, the following relations are generated 

M Tp y " = = ~ 2Tcusp (4,4) 

V 1 = =r “‘ ! ” (45) 

where T CMSp is the cusp anomalous dimensions of two light-like Wilson lines [37] . 

Figure (1) shows the path in (/x, u) space we have chosen to resum the logs. The hard 
part is run from the scales rrih down to the scale 1/bo ~ P. L> while the jet function is run in 
v space up to the scale rrih, where the scale bo is defined as be 1E / 2. In performing the hard 
resummation, n 2 are resummed by analytically continuing the matching scale from space-like 
to time-like kinematics relevant for Higgs production [39]. It has been shown that the inclusion 
of the 7r 2 in the resummation does improve the perturbative convergence of the hard matching 
expansion [19, 40], however this method only sums a subset of 7r 2 and one can not claim that 
this method leads to a systematic reduction in theory errors 5 . 

At NNLL, the systematics require that we keep the cusp anomalous dimensions at three 
loops, the non-cusp at two loops and matching at both the high and low scales at order a. All 

5 There is no scaling rule for these terms 
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Figure 1. The path in (/x, ix) space used to resum the logs. 


of the necessary anomalous dimensions have been calculated previously in the literature. The 
hard anomalous dimensions can be extracted from [43]. We utilize the one loop matching for 
the jet and soft function from ([15]). To run the jet function in u we need only the two loop 
cusp piece since log in jet rapidity the anomalous dimensionvis not large. The coefficient of 
this log is fixed by the hard function anomalous dimensions. We also the need the two loop 
“non-cusp” piece in v which can be obtained from the two loop jet function. The results for 
the hard anomalous dimension are tabulated in [44]. 

The matching of the TMPDF onto the PDF may be written as 


f^pM-O = Yl~ z j z 


+ 


( P±P± | 9± 
V p'i 2 


I± 2 g/k(z/z',$_)}fk/p(z') (4.6) 


IP-lI 


where the sum is on species of partons, and the parton PDF’s are defined as 

fg/p{z) = -zn ■ PnO{z) g ±/lu (p n \ [B^ ± (0)(5(z n ■ p n - P)B% L (0)] \p n ). 


fq/p(z) = -zn- PnO(z) (Pn 


e p (0 )W n U(zn ■ p n - V)w£&J0) 


Pn)- 


(4.7) 


We adopt the mostly minus metric such that conventions that p c [p^_g_Lap = — p\ and make use 
of the’t Hooft-Veltmann scheme, so the external transverse momenta remains in 2 dimensions, 
as do the external polarizations on the operator (the free Lorentz induces). The scheme choice 
is advantageous, as it allows one to renormalize the operator directly in p± space. At tree 











level the TMDPDF matches on to the PDF as follows: 


f (0)a/3 
•hi - 9/9 


(z,p±) = 5(1 


- z)5 { ' 2 \p 1 _) 



fg/p{ z )- 


(4.8) 

(4.9) 


The one loop renormalized TMPDF, necessary to extract the order a matching coefficients 

(-ffig/g)’ is S iven b y 


1 lAglg(zld,A) = ~~~^0 ^ “ ln “ z ) + Pgg*( z ) 

r ( P± 

AoUb — 


/ -2 \ n Ot s C A 1- z 


I±2g/g( z / z ,P±) = - 2 


7r 2 


/ P 1 P 1 

’ p ' v p± 2 J 


I±lg/q( z / z ,P±) 


ql s Cf 


7 r 


£o(/u,^)Pg 9 (^ + 5 (2) (p±)^(^) 


(4.10) 


where we have written the expression in terms of plus distribution C n = 
the properties of this distribution are collected in the appendix of [15]. 

The resummation is carried out in impact parameter space and when Fourier transforming 
back to p± space, the integrand diverges along the path of integration due to the Landau pole 
at 1/6 ~ A qcd- However looking at the form of the cross section, the integrand over 6 
contains the Bessel function of first kind with argument 6pj_ and an exponent which contains 
powers of ~ log(Mhb). The combination of these two factors quickly damps out the integrand 
outside the region l/p± > b > l/m,h. Thus to avoid hitting the pole, we can simply put a 
sharp cut off for 6 at a value sufficiently higher than the largest value of l/p± ~ 0.5GeH _1 
but smaller than 1/A qcd- A convenient value is found to be 6 2GeV 1 . It has been 
numerically ascertained that the variation of this cut-off between the values of 1.5 — 3 GeV -1 
produces an error only in the fourth significant digit of the integral. 



5 The Fixed Order Cross Section and Power Corrections 

When we match the full theory onto SCET we drop terms which are beyond leading order 
in A = p±/Mh. However, a comparison of the full theory fixed order cross section([10]) with 
the effective theory calculation reveals that in order to achieve ten percent accuracy, we need 
to include these higher order terms in A beyond p± > 30 GeV. In principle we could achieve 
this by keeping higher order terms in the factorization theorem. However, this would only be 
necessary if we wished to resum the logs associated with these power corrections. Given that 
these power corrections are only relevant in regions of larger p±, these logs are not numerically 
large and thus represent a small corrections which need not be resummed. At large p± we 
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should turn off the resummation and use the fixed order NNLO result. Thus we subtract 
out the logs from the NNLO correction so that we do not double count, and turn off the 
resummation using scale profiles which we now consider. 


6 Profiles 


The SCET formalism we have used retains the leading order operators in the power counting 
parameter A. This restricts the validity of the resummed cross-section to the regime where 
p _l << My. At larger values of of p±/Mf l the resummation becomes superfluous. Moreover, 
the power corrections from subleading operators become increasingly relevant and to maintain 
accuracy it is vital to include these power corrections, as discussed above. So it is desirable 
to smoothly switch over from the resummed result to the full theory fixed order cross section 
when the singular and non-singular terms of the cross-section are comparable. Notice that 
formally, it is not necessary to turn off the resummation. As long as one makes the proper 
subtraction to avoid double counting, then the resummed terms that are kept should be sub¬ 
leading in the logarithmic power counting, since we are in the transition region. However, as 
discussed in the introduction, the fixed order result involves the cancellation between singular 
and non singular terms. If the scale of the log is not chosen appropriately (of order p±) then 
this cancellation is not effective and the cross section is over estimated beyond the canonical 
scale variation errors 6 . This will be demonstrated below. 

The shutting down of the resununation is achieved by varying the low matching scale 
from its value 1/bo i n the low pj_ region to My in the high p± region. This has the effect of 
turning off the resummed exponent. To do this, we introduce profiles in p, v following the 
work done in [7, 8]. We use three typical profiles for varying the renormalization scale from 
l/6o t° Mh- Each profile is chosen as a linear combination of hyperbolic tangent functions 
which smoothly transitions from the resummation region to the fixed order one. The general 
equation for each profile P can be written as 


P(s,L,H,P,t ) 


L 

2 


^1 — tanh 


s 


AP 

t 




^1+tanh 


s 




( 6 . 1 ) 


where s determines the rate of transition, L = 2 e~ lE /b is the initial value and H = My is the 
final value, t is the value about which the transition is centered. Fig. (2) shows three profiles 
with s = 2 and t = (35,45, 55) GeV. 


7 Error analysis 

In our formalism, we have performed an expansion of the cross-section in My/rrit , p±/My, 
A.qcd/p± and the strong coupling a s . Therefore, the dominant error in our cross-section is 
due to the first sub-leading term that we drop in each of these parameters. When expanding 

6 We thank Iain Stewart for discussions on this point. 








Figure 2. Profiles for turning off resummation 

in Mh/mt, we retain only the leading order operator. A comparison of the cross section in 
this limit with the exact mt dependent result [10, 11], reveals that this approximation works 
extremely well for p± < mt- For the range of p± discussed in this paper, the error due to these 
corrections is less than 1 %. For the parameter p_\_/Mh , again we retain only the leading 
order operator. This approximation works well for p± < 30 GeV as the leading correction 
scales as p\/M^. As discussed earlier, for larger values of p±, we include the full NNLO 
expression which includes all orders in p±/Mh- 

For the third parameter, Aqcd/p± one must consider the non-perturbative contributions. 
Unlike the case of a central veto [44] there is more then one relevant non-perturbative param¬ 
eter. The reason is that we are working in SCETII where the collinear and soft function both 
have the same invariant mass scale. As such, there will be non-perturbative contributions 
coming from both of these sectors. As p± approaches A qcd-, we need to systematically include 
higher dimensional operators in both the soft and the beam sectors. There are multiple non- 
perturbative functions which contribute in this regime. In principle these matrix elements 
can be fit to the data. Work in this direction has been performed for vector boson production 
where ansatzes [47] for these matrix elements have been fit to the data [45]. However, since 
we are interested in a gluon initiated process we can not hope to use any extraction from the 
quark initiated vector boson production process. Therefore we make no attempt to model 
these corrections and simply include a rough estimate of the errors due to these terms by 
assuming an error which scales as A qqd/p±- 

Finally we consider the errors due to higher order perturbative corrections. The accepted 
methodology for estimation of these errors is to vary the scale p by factors of two. The idea 
is that such variations estimate the size of the constant terms which are not captured by 
the renormalization group, as well as the sensitivity to the perturbative trunctation of the 
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Figure 3. low px scale variation 

anomalous dimensions. Of course this is just a rough guess which falls into the rubric of “you 
do the best you can,” though it is important to try to gauge all sources of large logarithms 
in these variations. 

When working with the RRG we generalize this technique since we now have two distinct 
types of factorizations, virtuality and rapidity, and within the RRG there are therefore two 
distinct exponentiations. In both cases there is a choice as to what sub-leading terms should 
go into the exponent. Varying the scales p, v corresponds to varying the size of these sub¬ 
leading pieces. Given that virtuality and rapidity factorization are independent mechanisms 
we should vary both of these scales in order to estimate the errors in the choice of exponents. 

This variation is accomplished by adjusting the parameter L in each profile which sets the 
scale for the low scale matching of the jet and soft function that dominates the perturbative 
errors. In keeping with the canonical recipe adopted by the community we vary L by a factor 
of two about its central value(l/6o). The corresponding modification in the profiles for either 
p or v is shown in the Fig. 3. When we vary these scales we must keep in mind the restriction 
that the argument of each large log is varied at most by a factor of two. This gives us a 
constraint 0.5 < p/v < 2 since there are logs of the form In {p/v) in the exponent. In all, 
this provides us with six different possibilities for each of the three profiles (Fig. 2). We 
do a similar analysis for the end point region of each profile by varying the scale H (Fig. 4). 
Since this region of transverse momentum is dominated by the fixed order cross section, the 
variation in H amounts to a fixed order scale variation. We choose to keep the largest error 
band generated due to these variations. The plots which combine all the effects above, are 
shown in Fig. 5 and Fig.6 for 5 = 8 and 14 TeV. We use the MSTW 2008 pdf at NNLO [51] 
and the value a s (M z ) = 0.1184. 


10 








Figure 4. High pr scale variation 



Figure 5. Cross section at 8 TeV 

We also consider how the cross section is typically affected by using other pdf sets. In 
particular we evaluate the cross section at 13 TeV with the CTEQ5 pdf. The deviation form 
the MSTW2008 result is of the order of a few percent as shown in Fig. 7. We include these 
errors as well in the final plot for the cross section at 13 TeV Fig. 8. 
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8 Comparison to Previous Results 


Let us now compare our results with the previous results in the literature which include the 
NNLL resummation. We will compare only with those papers which show plots which include 
all of the corrections considered here, i.e. [21, 52], We agree formally with the results of these 
papers. Of course formal agreement is not necessarily the relevant issue. Formal agreement 
means that the coefficients of the logs in the resummed result (as well as the non-singular 
pieces) agree. However, it does not mean that the argument of the logs are the same. Indeed 
for sufficiently different arguments the results can differ by an amount which is well outside 
any “reasonable” theoretical error estimation. 

In the work [21] the authors calculate the spectrum up to 60 GeV and match onto the NLO 
result. Terms of order a 2 are not included which is consistent as long as one presumes that 
^ln(60 2 /mh 2 ) ~ 1. That is, as long as the logs still dominate all the way up to p± = 60 GeV , 
then keeping the terms of order a 2 would not improve the accuracy of the calculation. The 
results in [21] dont track the rapidity scale. Meaning that an implicit fixed matching scale 
has been chosen. In contrast to [21], the NNLL error band in our analysis lies within the 
error band of NLL. The reason is mainly two fold. They have no rapidity scale to vary and 
they have no errors due to non-perturbative corrections. This elucidates the importance of 
including the scale v in error analysis for a better estimation of accuracy at each order in 
resummation. Further, the authors claim the existence of a non-perturbatively generated 
hard scale on the order of 8 GeV, which they state cuts off the non-perturbative corrections 
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Figure 6. Cross section at 14 TeV 
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Figure 7. Comparison between pdf sets MSTW2008 and CTEQ05. The y axis plots the fraction of 
the difference between the central values of the cross section using the two sets 
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Figure 8. Cross section at 13 TeV 
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Figure 9. Comparison with Becher et. al. 



Figure 10. Comparison with CSS formalism 
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and therefore there results are valid down to vanishing p±. The underlying justification is 
based on the paper [50] where it was argued that since in QED the probability to emit a 
photon with p± ~ 0 is vanishes faster then any value of Q 2 , the only contribution to the rate 
must come from two photons whose total p± approximately vanishes but whose individual 
p _l is large. The authors of [50] then conjecture that the same reasoning should hold in the 
non-Abelian case. A recent study of the non-peturbative contributions to the TMPDF can 
be found in [55] and we invite the reader to consult this paper for further information on the 
validity of these claims. Fig. 9 shows a comparison with their central value at NNLL+NLO. 
For larger values of p±, the central value curve is below our estimates since the matching is 
done with fixed order NLO. 

Grazzini et. al. [52] use the CSS formalism and plot the result out to 120 GeV, matching 
onto the fixed order result at NNLO. They do include an additional scale they call the resum¬ 
mation scale Q, in addition to the factorization scale of the PDFs and the renormalization 
scale of a. The variation of Q variation mixes what we would call p and v radiation. This is 
discussed in detail in the appendix. To turn off the resummation in the exponent they make 
the replacement 

Log(Q 2 b 2 /b 0 ) ->• Log(Q 2 b 2 /b 0 + 1) (8.1) 

this has the effect of killing the resununation at asymptotically small b and insuring that the 
total inclusive cross section in reproduced. However, in the range Mh/p_\_ ~ 1 the resummation 
is not shut off rapidly enough. Indeed, our results disagree with [52] for the central values 
by an amount which is larger then the theory error bars in the tail region as can be seen in 
Fig. 10. The reasons for this overshoot, as discussed above, is that if one does not use a 
profile to shut off the resummation beyond the transition region, then the scale of the log in 
the large p± region is such that the singular and non-singular terms in the fixed order cross 
section no longer cancel as they do in the fixed order result. This effect is a recurring theme 
in the literature [44, 48, 49]. The authors of [52] also include the two loop matching to the 
hard function and show that its effects are at the one percent level. Moreover, without the 
other pieces of the calculation of this order, 4-loop cusp, 3 loop non-cusp, two loop soft and 
collinear matching, this inclusion does not increase the accuracy of the calculation. Regarding 
the scale variation, the authors vary three distinct scales as was performed in this paper. The 
variations mix rapidity and renormalization group scale dependence (see appendix) and one 
could argue that this is perhaps not as clean as varying p and v independently. But given 
that the scale variation technique is a blunt instrument to estimate perturbative uncertainty, 
it is not clear that distinguishing logs is quantitatively relevant, as long as all resummed logs 
are probed in the variation. The uncertainty bands in [52] are very close to those presented 
here. In the appendix we show how the work of Grazzini et. al. can be parsed in terms of 
the RRG. 

Finally, a recent paper [54] working within the TMDPDF formalism of [24, 25, 27] pub¬ 
lished a cross-section for the Higgs transverse momentum distribution at NNLL, with no fixed 
order matching. The error bands are very small, with no overlap between the NLL and the 
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NNLL resummation, which the authors noted as a sign that they had under-estimated the 
perturbative uncertainty. The bands were generated by varying the resummation scales gen¬ 
erated by the /r evolution of the TMDPDFs, after the resummation of the rapidity logs, which 
followed the procedure of [56]. The rapidity resummation scales were taken as fixed quanti¬ 
ties. Thus no attempt was made to directly gauge the perturbative uncertainty in the rapidity 
resummation, that is, the perturbative uncertainty in the Collins-Soper kernel, whose expo¬ 
nentiation also resums the rapidity logarithms. Interestingly, variation of the non-perturbative 
models for the Collins-Soper kernel produced a much larger impact on the cross-section, see 
Fig. 4 of [54]. Since these non-perturbative models are explicitly exponentiated with the 
perturbative Collins-Soper kernel, in our view, this variation of the non-perturbative model 
parameters then may not be an estimate of the non-perturbative physics at all, but a more 
indirect gauge of the perturbative truncation of the Collins-Soper kernel and the rapidity 
resummation. Separating out the effects of the perturbative uncertainty in the Collins-Soper 
kernel, which has not been studied in the literature, will prove important in isolating the 
non-perturbative physics of the TMDPDFs. 

9 Conclusion 

The next LHC runs promises to shed more light on the nature of the Higgs boson. Here we 
have attempted to give the most accurate result possible for the transverse Higgs spectrum 
using the currently available theoretical calculations for matching coefficients, anomalous 
dimensions as well as the fixed order cross section. We resum the large logs in the low 
p j_ region to NNLL and match this result to the fixed order cross section for high values 
of transverse momentum. To maintain accuracy over the complete range of the transverse 
momentum, we have utilized profiles in both the virtuality (fi) and rapidity {y) factorization 
scales which smoothly turns off the resummation and matches onto the NNLO fixed order 
cross section at large values of p±. This procedure prevents the problem of enhancing the 
fixed order result. We discuss the sources of error and ways to effectively estimate them using 
possible variations in the rapidity and renormalization scales for the profiles introduced. While 
we get a good agreement with previous results in the low p± region, we get a lower estimate 
for the cross section in the transition and tail region. 

The code which generated the plots in this paper is available upon request. 
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A RRG resummation 


In this appendix, we review the RRG formalism for the purpose of explicitly comparing to the 
method used in [31, 32, 53]. After integrating over the rapidity of the Higgs, the differential 
cross section is given by 

-7-2 = a o f —— f d 2 be l ^S(mh- Sxix 2 )H(m h , p) 
dpf J xi x 2 J 

fn±(^lX fn±(x2,b, IJ,V^s(b, , (A.l) 

where we have written out the explicit momentum fraction dependence. Since we are focusing 
on the resummation, we will often suppress the other arguments of the renormalized functions. 
The soft and beam functions have the /z-anomalous dimension: 


M-y-ln S(u,p) = r R [a s (^)]ln[-l +7 s (a s (p)), 
dp Lz/J 


d 


1, 




lm h J 


+ 7 /± (a s (//» 


(A.2) 

(A.3) 


Using the fact that the RG and RRG variations of the amplitude must commute, we can 
write the rapidity anomalous dimension as integral over the p anomalous dimensions [14, 15] 
and a constant piece 


d f 1 

p—lnS(u, p) = / 

du v J(^y Q 


-1 +7i?fa s (( 6 ^) X ) 




dq 


1 


,be^ E s -1 


J T ^s(q)}-^ R [as(( — r ) 


(A.4) 

(A.5) 


Writing things in this way allows for there to be large p dependent logs which may arise 
depending upon the choice of p. As noted in Sec. 4, Tr is related to the cusp anomalous 


can be 
-1 


dimension. The non-cusp piece of the rapidity anomalous dimension 7 r 07 (( 7 ^) ) 

extracted from the calculation of the soft function at the renormalization point p = ( 77 ) 
This is a choice of scheme which defines the scale in the Log which multiplies the cusp 
anomalous dimension. 

We run the beam/soft functions in v from v s and m/,, to v respectively 


S (b, p, = U^ RG (b, p, S (b, p, z/ s ) , (A. 6 ) 

fn±(x,b,p,i^j = U RRG (b, p, ^-) f n ±(x,b, p,mhj . (A.7) 

Note that we have explicitly chosen vs = rrih in the TMDPDF. The full RRG factor has the 
rapidity factorization scale v cancel in the exponent, and gives: 


tt RRG 

u sf 2 



—) = Exp 
Z7 J 


Inf?* 

V v 



(1™)-. 7 r HKW] + to 


a. 




(A. 8 ) 
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We also need the /j evolution of the TMDPDFs at v = mh and the soft function at v = v s . 
With these scale choices, the TMDPDF has no double log-/x evolution, so one obtains: 

f±(x, b, n, v = m h ) = Exp 


/ 77/1 [«*(?)] 

' Hi " 


= m h ) 


S(b, H, V = i/ s ) = Exp 


r 11 d q< 


Li q A 

l V q J 


S{b,Hi,v = v s ). (A.9) 


Combining these pieces together, we achieve the full resummation kernel for the cross-section: 
/nJ_0, v)fn±-{v, v)S{h, v) = R[h, Hi, v s ,m h , bjf n ±(fj.i, m h )fn±(m, m h )S(m, v s ) (A.10) 

+'ys[a s (q)] + 27/±[a s (g)]} 


) = Exp 

[f In 

[m h l 

/ 

[ Jim c l 1 

1 q \ 


+ In 


m h 


v. J 


f dq n r / m r w6e 7s \_i,i 

r « [ “' ( « )1 + ™[“*<(—) >]} 


(A-ll) 


From the form of (A. 11), we can directly see the variation of the independent resummation 
scales fi->Hii u s probe different exponentiated logs, thus giving estimates of the subleading 
terms in the perturbative expansions. The double logarithmic terms associated with fii cancel 
manifestly in the exponent', and thus the Hi variation is estimating the subleading terms 
associated with 7 s j ± - We can eliminate large logs of the impact parameter from the beam 
and soft function matrix elements by choosing the scales Hi ■, v s appropriately. Four potential 
schemes for canonical scale choices are: 

A : Hi = u s = 

B : Hi = PT, 1 

C : Hi = = PT 

D : Hi = v s = Pt 

As noted in [57], making scale choices in either momentum space or conjugate space can have 
a sizable numerical impact on the cross-section, as well as how accurately the resummation 
captures the higher order logs claimed by the resummation. Scale setting in momentum space, 
residual logs can exist in the plus distributions of the low scale theory that become apparent 
only after integrating the resummation kernel against the matrix elements, in this case, the 

'Though depending on how one handles the integration over the running coupling, this cancellation could 
be incomplete. An incomplete cancellation could be used to give another handle on estimating the uncertainty 
in the double logarithmic terms. 
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TMDPDF and the soft function. Scale setting in conjugate space does not suffer from this 
ambiguity. Specifically, for the canonical choice u s = /i, = ( he 2 K ) , R becomes: 


R[n,m h ,b) = Exp 




|W) 1 1 

+ In 


mh 
q J 


r fl[«s(g)] + 7s [«*(<?)] + 2 7 /j.N 


m h be lE 7 


(A.12) 


A.l Comparision to Grazzini et al. 

We now write out the resummation formalism used in [52], as derived in [31, 32, 53], sup¬ 
pressing the flavor sum in the PDFs: 


da 

% 


dli dl 2 f ( b,m h , XlX2 , fi R ,/J, F 


dpf J Xl X 2 J 47T 


(A.13) 


Then one switches to moment space by integrating A = X 1 X 2 for W: 


fi 


W N \b,m h ,p R ,p F j=J dz z N b,m h , ^,/j, R ,p, F 


1 


(A.14) 


The authors write this as 


Wat | 'b , m h , p R , p F ") = R N (rn h , p R , p F ] Exp 


G n [ b, m h , fj, R , Q 


(A.15) 


Note that the factorization scale from the PDFs is in the function R. Then the form of the 
resummed exponent is: 


Gn\ b, m h , p R , Q = - 


r d 4\ A 

lb 2e 2 TE <T 


a s (q) 


ml 


In [-±)+B N 


a s{q) 


(A.16) 


Note that Bn has PDF running in it since it depends on the moment N. The function Rn 
is claimed to have no large logarithms, and can be perturbatively calculated. To one loop: 

U N (m H , hf; qJ = 1 + ^ + 2Ctf - pfSoin + 2 lN t F - (*^ ( % + £ (1) + 2 7 ^)^ 


(A.17) 


Where: 


£r = ln^ 

A4 


l F = In ^ 
d'F 


£q = ln^f (A.18) 
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And: 


A (1) = 

= - l -(llC A -2N f ) 

A> - '* 


dzz N ~ l P {l \z) 


(A.19) 
(A.20) 

(A-21) 


Cat is related to the order e pieces of the DGLAP splitting kernels. H W is the hard matching, 
but also contain terms in the 5(1 — z) pieces of the TMDPDF. Note that hr ~ hf ~ Q ~ m/i, 
so that the PDFs are run from A qcd to the high scale, implicitly. 

By choosing the appropriate scales, we can connect the resummed formula of (A. 16) to 
the resummation formulas of (A. 11) and (A. 12). We must also apply the resummation scale 
prescription of [53], where all large logs are split: 

\m h - 


lnK(V) ^ 


In 


lq J 


+ In 


Qi^r 1 


(A.22) 


and any log of In 


m h 

Q 


is expanded out of the exponent and included in the hard function of 


(A. 15). Setting fi = Q in (A. 12), this gives: 


R(Q,mh,b ) = Exp 


Q 

m k 


lijln 

Q 1 


m h 


q J 


r^Mi)] + 7« [««(?)] + 2 vA a M)\ 


+ In 


VQ\ 


7 Rias^y 1 )] 


r (( helE \- 1 
7 r[<xs{{-z-) 


(A.23) 


choosing Q = fi this result is identical to (A. 12). 

The hard function of Eq. (A. 15) is explicitly given to two loops in [53] , and one can 
verify the factor In ^ 7 _r[ a s((^yr^) X )] i n (A.23) is expanded out of the exponent and then 
included in the B ^ term of the hard matching. Indeed, from the form of the hard function 
given in [53], one can see the resummation scale Q probes both rapidity logs and RG logs 
since B ^ includes both types of anomalous dimensions. This explains the similar sized 
perturbative error estimation of this paper, and that of [52]. 
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